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Abstract. In Newtonian gravity, a stationary axisymmetric system admits a third, 
Carter-like constant of motion if its mass multipole moments are related to each other 
in exactly the same manner as for the Kerr black-hole spacctime. The Newtonian 
source with this property consists of two point masses at rest a fixed distance apart. 
The intcgrability of motion about this source was first studied in the 1760s by Eulcr. 
We show that the general relativistic analogue of the Euler problem, the Bach-Weyl 
solution, does not admit a Carter-like constant of motion, first, by showing that it does 
not possess a non-trivial Killing tensor, and secondly, by showing that the existence 
of a Carter-like constant for the two-center problem fails at the first post-Newtonian 
order. 
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1. Introduction and summary 

Every physics undergraduate knows that solving dynamical problems in classical 
mechanics is simplified if one can find constants of the motion. Given enough suitable 
constants of motion, the problem becomes completely integrable, solvable in terms 
of integrals. The simplest example is a static, spherically symmetric potential. In 
Newtonian gravitation, the conservation of energy, total angular momentum and one of 
its components guarantees that the motion for this problem is completely integrable. 
Likewise, in general relativity, given the same symmetries, the conservation of the rest 
mass (or the unit norm of the four velocity) in addition to the other three constants of 
motion leads to an integrable system. But for problems with less symmetry, completely 
integrable systems are rare. 

Systems that are stationary and axisymmetric possess conserved energy and angular 
momentum about the symmetry axis (and conserved rest mass in the relativistic case). 
But under what conditions do they possess an additional constant of the motion, 
sufficient to make them integrable? Recently, one of us [1] pointed out that a stationary 
axisymmetric system in Newtonian gravitational theory could possess such a third 
constant of the motion, analogous to the Carter constant of the Kerr geometry, if the 
multipole moments Qi of the source obey the same constraints as those satisfied by the 
mass moments of the Kerr spacetime, namely 

Q 2e = ma 21 , Q 2e+1 = , (1) 

where I = 0, 1, 2, . . .. Unlike the Kerr case, where a is the angular momentum per 
unit mass, here a = (Q2/m) 1 ^ 2 , where Q2 is the quadrupole moment and m is the total 
mass of the system. Nevertheless the similarity of Eq. ([I]) to the Kerr sequence of mass 
moments is striking. 

In the prolate case, Q 2 > 0, the Newtonian source that has this property 
consists of two bodies each of mass m/2 held at fixed locations ±a on the z-axis. 
Subsequent to the publication of Ref. p], its author learned that this problem was 
already studied in the 18th century by none other than Leonhard Euler, and is now 
known to classical dynamicists as the "Euler problem" , a non-spherically symmetric, yet 
completely integrable system (for a thorough review of the history and mathematics of 
the Euler problem and related integrable systems, see [2]). Other authors have similarly 
rediscovered the unique properties of the Euler problem in other contexts, including 
the quantum mechanics of the singly ionized hydrogen molecule [3] (the underlying 
electrodynamics is mathematically equivalent to the gravitational Euler problem) and 
stellar dynamics [I]. 

In the oblate case Q2 < 0, a is imaginary, and the source consists formally of two 
equal point masses at ±ia on the imaginary z-axis, although the gravitational potential 
is still perfectly real. 

It turns out that the solution of Ref. [1] can be generalized to the case where the 
dipole moment Q\ is non-zero; here the sources in the prolate case are two unequal 
masses held at ±a on the z-axis, with m + + m„ = m. In the oblate case, the sources are 
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two complex conjugate masses on the imaginary z-axis. Some of these generalizations 
had previously been studied by celestial mechanician John P. Vinti [5] , who noticed that 
if the parameter a of the Euler problem was adjusted to match the Earth's Q2, then the 
value of Qi implied by Eq. (pjj closely approximated that of the Earth (though higher 
moments did not), so that the integral solutions of the Euler problem gave reasonable 
approximations for orbits of satellites around the Earth [2] . 

These results in Newtonian gravity raise the obvious question: is there an analogous 
third integral of the motion for the general relativistic version of the Euler problem, 
namely the Bach-Weyl solution [6J [T] ? This is an exact, static, axisymmetric, vacuum 
solution of Einstein's equations, corresponding to two point masses held a fixed distance 
apart on the z-axis. It is asymptotically flat, and is regular everywhere except along a 
line joining the two masses; on this line there is a pressure singularity representing a 
"strut" needed to hold the masses apart. 

In this paper, we show that in fact the Bach-Weyl solution does not admit a Carter- 
like constant analogous to that of the Newtonian Euler problem. We do this in two ways. 
First we show that the Bach-Weyl spacetime does not admit a non-trivial second-rank 
Killing tensor £ a p. Had it done so, then by virtue of the Killing equations C(afi;-y) = 
and the geodesic equation p^p a .p = 0, there would have existed a conserved quantity 
given by C = £, a /3P a P^ ' ■ It is important to recall that the metric itself and any product of 
Killing vectors are also Killing tensors, but they are trivial in that they do not generate 
new constants of the motion. 

Second, we show that, in the post-Newtonian limit of the Bach-Weyl geometry, a 
Carter-like constant analogous to that which exists at Newtonian order cannot be found. 
Whereas in the Newtonian Euler problem, the Carter-like constant is given by 

C = h 2 + a 2 v 2 z - 2azU* . (2) 

where h = x x v is the angular momentum per unit mass, and 

^^±_^h (3) 

where r± = |a? =f ae z \ and m± are the two masses (in general they do not need to be 
equal), in the post- Newtonian Euler problem, we show that the closest one can get to a 
conserved quantity is the expression 

1 d 
2~dt 

where U = Gm + /r + + Gm_/r„. The term on the right-hand side cannot be expressed 
as a total time derivative of any function of x or v, and so a Carter- like constant of the 
chosen form does not exist at post-Newtonian order. 

The remainder of this paper gives the details underlying these conclusions. In Sec. 
II we discuss the Bach-Weyl solution and its symmetries. In Sec. Ill we derive the 
general solution of the Killing tensor equations, and show that the result is a linear 
combination of trivial Killing tensors. Section IV takes the post-Newtonian limit of the 
Bach-Weyl solution and shows by construction that a post-Newtonian generalization of 



8U\ J v 2 AU 



6 TT* dU ,A^ 
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the Euler-problem's Carter-like constant does not exist. Concluding remarks are made 
in Sec. V. In an Appendix, we give a more detailed argument supporting a key step in 
the analysis. 



2. The Bach-Weyl solution 

The Bach-Weyl metric [6], [7J, [8] is given in cylindrical p — z coordinates by 

ds 2 = -e 2X dt 2 + e- 2X [e 2u (dp 2 + dz 2 ) + p 2 dtf} , (5) 



A(p,z)=-^±--, (6) 
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where 



m 2 _\ m + m_ ( p 2 + z 2 - a 2 \ 

where r\ = p 2 + (z =)= a) 2 . Here and henceforth, we use units in which G = c = 1. 
The total Kepler mass of the system is given by m + + m_. The physical system can be 
viewed as consisting of two point masses at z — —a and z = +a; the metric is regular 
everywhere except at the two points, but the Einstein tensor (and thus the stress energy 
tensor) diverges along a line joining the two points, representing the "strut" required to 
hold the two masses a fixed distance apart. 

The metric admits the timelike and azimuthal Killing vectors, = d/dt and 
£W = d/d(f>. If it also admits a symmetric, second-rank Killing tensor £ a p, then 
satisfies the equations [8] 

Ka/3-y = ^7/3;a + C«7;/3 + Ca/3;7 = , (8) 

where Greek indices run over the values ... 3, semicolons in subscripts denote covariant 
derivatives, and the notation K a ^ labels the partial differential equations. The metric 
g a p and any symmetrized product C \ a C fi) °f Killing vectors automatically satisfy 
the Killing tensor equations, and thus our goal is to find non-trivial solutions of Eq. (IS]). 
Only a non-trivial Killing tensor will generate a new constant £, a pV a V °f the motion. 



3. Searching for a non-trivial Killing tensor 

We define a template for the symmetric Killing tensor as 



( A 
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Q \ 
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(9) 



where the 10 functions A,...,Q are a priori functions of t, p, z and 0. Notice that, 
notwithstanding the fact that the metric is stationary and axisymmetric, it is incorrect 
to assume a priori that the Killing tensor is independent of t and 0; in spherical 
symmetry for example, the trivial Killing tensors constructed from products of the 
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two non-axial rotational Killing vectors depend explicitly on 0. Using Eqs. fl5]) and 
in Eq. (jSJ), we obtain a system of 20 coupled partial differential equations. Following a 
method suggested by Garfinkle and Glass [9] we divide the equations into three groups 
corresponding to the (p, z) block (the functions D, C and E), the (p, z, 0) block (H, 
J and B) and the remaining (t, p, z, 0) block (A, E, G and Q). We define a = A — v. 
(a) (p, z) block: 

^.o-f^^pa-^g^, do, 

Jf„:0 = ^+2C^-2B^. (13) 



<9z <9z 9p 



(b) (p, 2, 0) block: 



dH / dX da\ da H ldD 
AW:0 =- + ff (2- + -j-J--2- + --, (14) 



(c) (p, 2, 0, i) block: 



dz dp 
dH dJ 

d\ da 

dp dp 

9 „ fdB d\ B dH 
^ :0 = e -(- + 4 B --4- + 2- 





da \ 






«9A 




<9z 






+ 90 





2l V<9B , „ D dA , <9J\ ^ 2 d\ 



(16) 



z 9A 1\ d\ i , 

2D " (>")- 2£ " sI- (17) 



^ : = e^- + 4 B - + 2-j-2<V- 

-2E^f^-iV (18) 



V«9p p 

^ :0 = ff (|-I) +J |-I^. (10, 



^ :0 = ^- F U + ^J- G U-^J + 2^' (20) 

K ^-- = ^- G {Yz + d-z)- F {dp--d-p) + 2m^ (21) 

K zt : = ^ + ^-2F^-2G^ + ^, (22) 
dA t A dX n 4A _ 2l/ f^dX , ^9A\ , O 0F 



*- : = a? " ~ 2e {% + %J + 2 aT • < 23) 
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AVi:0 =2^-2^^- + J-j + -e-, (28) 
9A 1\. , „3A _ 2 3Q 2 „ 1 _ 2 9B 2 „ 



: = F [ir P —,) + G w, - » ir/ v - tT sT ■ (29) 

First, we focus on the p, z block, Eqs. ( fTOi) - ( fT3l) . After changing variables to 

C = W + F\/ 2 , D = WV -YV 2 , E = X\/ 2 , (30) 

where V = e _2a , and W, X and F are functions of t, p, z and 0, we obtain 

K pzz : = W,„ + FpF - FKp + 2X,y + , (31) 
K zpp :0 = W, Z - Y, Z V + YV Z + 2X, p V + XK P , (32) 
K ppp :0 = W tP - Y p V - YV P + XV z , (33) 

k zzz : o = w% + r 2 v + *%, + xv p , (34) 

where commas in subscripts denote partial derivatives. Subtracting Eq. ( 133]) from Eq. 
fl3D and Eq. ([32} from Eq. ([34]) yields 

r p + x, = 0, r z -x p = o, (35) 

whose general solution is 

Y + iX = F(£,</>,t), (36) 

where J 7 is an arbitrary analytic function of the complex variable ^ = z + ip = re l9 , 
and an arbitrary function of the variables <fi and t. Notice that J 7 is proportional to the 
variable t defined by Brink [TO] . 

Taking the first derivative of Eqs. (]3Tj) and ( 132]) with respect to z and p respectively, 
subtracting, and using the fact that d/dp = i(d/d^ — d/d£) and d/dz = d/dS, + d/d£, 
we obtain 

0, *) + 3K?(^, 0, *) + 2K«(^, 0-F(£, 0, *) = cc. , (37) 

where cc. denotes the complex conjugate of the left-hand-side. This equation can also 
be written in the form 

'CP*), — . 

This is equivalent to Brink's Eq. (16). 
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Because V = e~ 2a( - p ' z ^ is a complicated, non-analytic function of p and z (or of £ 
and £), we can conclude that Eq. fl37|) or f )38|) can be satisfied for arbitrary p and z if 
and only if 

^(£,<M) = 0, (39) 

which immediately gives 

X = Y = . (40) 

A more detailed argument supporting this claim is presented in the Appendix. From Eqs. 
(J3T|) and f[3"2"j) we then conclude that is purely a function of 4> and t i.e. W = Ci(0, t), 
and thus that 

C = D = Ci((f), t)e~ 2a , £7 = 0. (41) 

Armed with the solution for C, D, and E, we now consider Eqs. (fT4|) - (TT9|) involving 
the functions if, J and B. After changing variables to 

J = jp 2 e- 2 ^\H = hp 2 e- 2 ^\B = C l {<P,t)p 2 e~ 2X + a iP i e- i \ (42) 

where j, h and C 3 are functions of t, p, z and 0, these equations become 








= h tP - ha, p - ja s + ^p 2 e 2X C w , 


(43) 







I 

= j, z - ha ;P - ja jZ + -p~ 2 e 2X C li(j} , 


(44) 


K^zp : 





= h z + j tP , 


(45) 


Kp^j, 





= 2h^ + p 2 C 3:P , 


(46) 







= 2j^ + p 2 C 3 , z , 


(47) 







= h(\ p - p- 1 ) + j\ z - -C 3 , - -p- 2 e 2A C 1)0 . 


(48) 



Combining Eq. (J45J) with the difference between Eqs. (143]) and fl44j) we obtain the system 
j, P + h z = 0, j,z ~h p = 0, (49) 

which implies 

j + ih = g(£,<l>,t), (50) 

where Q(£,,(f),t) is an analytic function of £. Substituting Eq. (15T)]) into the sum of Eqs. 
( HUD and flU} gives 

[V{£, 0, t)U + (cc.) + C li0 (0, t)V(£, Op~ 2 e 2X ^) = . (51) 

Again, because A and v are complicated non-analytic functions of p and z, we conclude 
that 

£(£,<M) = 0, Ci,* = 0, (52) 

which immediately gives j = h = and thus J = H = 0. Accordingly, Eqs. fj4*6l) . f|4Tj) 
and (|4"8"|) . imply C 3jP = C 3jZ = = 0, and thus that 

B = C 1 (t)p 2 e~ 2X + C 3 (t)p 4 e- 4A . (53) 
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We can now consider the remaining 10 Killing tensor equations. After changing 
variables to 

A = _ Cl {t)e 2X + C 2 e iX , Q = C A p 2 , F = fe 2v , G = ge 2u , (54) 
where C 2 , C 4 , / and g are functions of t, p, z, and 0, Eqs. (120]) - ( 129]) take the form 



K pp t 





1 

= f, P ~ foe, p ~ gu, z + 2 C ' 1 .* e 


-2A 

5 


(55) 


K zzt 





= g, z - ga,z - fa p + 


-2A 

) 


(56) 


Kp Z t 





= g, P + f,z , 




(57) 


Kp tt . 





= 2f )t + e* x - 2v C 2tP , 




(58) 


K ztt : 





= 2g, t + e^C 2 , z , 




(59) 


Km : 





= f\ P + g\ z + \ {e- 2X C u - 


■ C2,t) , 


(60) 


K<f, p t 





= p 2 C 4 ,p + e 2 »l^ 




(61) 


K^ z t 





= P C 4)Z + e gj, 




(62) 


K^t 





= 2p 2 C 4 , t + e 4A C^, 




(63) 


K t <p4, 


: 


= / (A, P - p^ 1 ) + g\, z - 


-\{C u e- 2X + C z , P 2 e-^). 


(64) 



Combining Eq. ( IBTj) with the difference between Eqs. (155]) and (156]) we obtain the 
system 



g, P + f, z = o , #, 2 - 



o 



(65) 



which implies 

g + if = H&<f>,t), (66) 

where "H(£, 0, t) is an analytic function of £. Substituting Eq. (166]) into the sum of Eqs. 
( 155]) and ( 156]) gives 

0, t)U + (cc.) + CM t)Vfa i)e~ 2X ^ = . (67) 

As before, we conclude that 

H(£,<f>,t) = 0, C 1)t = 0, (68) 

which immediately gives / = g — and thus F = G = 0. Accordingly, Eqs. ( 158]) - (162j) 
imply C 2 , P = C 2jZ = C 2 j = 0, and C 4iP = C^ z = 0. Equations ( 163]) and ( 164]) then yield 

dC 4 (<M) 



<9C 4 (<M) 



24X dC 2 (<p) 



90 



.2 -4Ag^(g 



o, 

0. 



(69) 
(70) 



Eq. (169]) can be satisfied for all p, z, and t if and only if C 4jf = and C 2 ^ = 0. 
Similarly, Eq. (170]) can be satisfied for all p, z, and i if and only if C^ = and 
C%t = 0. The final result is that C±, C 2 , C3 and C 4 are arbitrary constants. 



Carter-like constants of motion 



9 



Combining these results, we find that the solution can be expressed in the form 

W = C l9aP + G&kf + C^Hf - 2C4$$ } , (71) 

which is a linear combination of trivial Killing tensors. We also note that this solution 
is obtained immediately using Maple. We conclude that the Bach-Weyl spacetime does 
not admit a non-trivial second-rank Killing tensor, and thus, unlike its Newtonian limit, 
does not admit an analogous Carter-like constant of the motion. 



4. The post-Newtonian limit 

Here we show by explicit construction that, in the first post-Newtonian approximation 
of the Bach-Weyl metric, a Carter-like constant does not exist. It suffices to note that 
the metric function u, being already quadratic in the masses, contributes to the motion 
only at 2PN order, and thus it can be set equal to zero in Eq. (JSJ). The result is a metric 
in isotropic coordinates which matches the standard form, say of the parametrized 
post-Newtonian (PPN) framework [11], with the GR values 7 = /3 — 1 of the PPN 
parameters. Since the metric is static, the post-Newtonian equations of motion for a 
test body take the form 

f?L = VU(1-4U + V 2 ) -Aw- VU, (72) 
at 2 

where U is the Newtonian gravitational potential given, in the Bach-Weyl case, by 
U = — A = m + /r + + m_/r_. Following the method used in the Newtonian Euler 
problem pQ, we calculate dh 2 /dt, where h = x x v: 

ldh 2 ( dv\ 

h • \ x x — 



2 dt V dt J 

= h- (x x W) (1 -4U + v 2 ) -Ah 2 ^, (73) 

where we use the fact that, for a stationary potential, v ■ 'VU = dU/dt. After some 
manipulation, it is straightforward to show that 

h-(xx VU) = ^(azU*) - a 2 v z V z U 
dt 

= -— (2azU* -a 2 v 2 ) 
2dt K zJ 

- a 2 v z (4<^ + AUV Z U - v 2 V z U^j , (74) 

where we used the PN equations of motion to go from V Z U in the first line to dv z /dt, 
and where U* = m + /r + — m_/r_. Considering just the Newtonian parts of Eqs. (|73|) 
and (I74p . we find the Newtonian Carter-like constant of the Euler problem, Eq. (j2]). 

Now including the post-Newtonian terms in Eqs. (1731 and (1741) . and using 
Newtonian equations of motion and the Newtonian Carter-like constant where necessary 
in those terms, we find 

~ [(h 2 + a 2 v 2 )(l + 8U) - 2azU*(l + v 2 + 4U)] = -QazU*^- . (75) 

& (JjL> (Alt' 
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The term on the right-hand side cannot be expressed as a total time derivative of any 
function of x or v, and so a Carter-like constant does not exist at post-Newtonian order. 

5. Discussion 

We have shown that a Carter-like constant of the motion analogous to the one that exists 
in the Newtonian Euler problem does not exist in the relativistic analogue given by the 
Bach-Weyl solution, both by showing the non-existence of a non-trivial second-rank 
Killing tensor and by showing the absence of such a constant at post-Newtonian order. 
However, this does not completely rule out the integrability of motion in the Bach- 
Weyl spacetime, because there could in principle exist a higher-order Killing tensor, or 
equivalently an additional constant of motion more complicated than the Carter-like 
constant we have considered. 

It is also worth pointing out that, because our solution did not depend crucially 
on the specific functional form of the function V, it applies to other solutions in the 
Bach-Weyl class, such as the two black hole solution (corresponding to a potential A 
generated by two collinear rods on the z-axis) discussed in [12]. In that solution, even 
the Newtonian limit fails to have a Carter-like constant. 
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Appendix: On the vanishing of the analytic function J 7 

In this Appendix we present an argument in support of the claim ( 1391) that the analytic 
function J 7 must vanish. Eq. fl37|) can also be written as a real differential equation in 
terms of Y and X, in the form 



We perform a Taylor expansion of the functions X, Y, and V about an arbitrary point 
(p , Zq), up to an order M, corresponding to the powers Ap q Az M ~ q , where q = . . . M, 
Ap = p — p , and Az = z — z$. Now because the function T = Y + iX is analytic, it 
can be expanded in the form 



2VY pz + 3{Y p V z + Y Z V P ) + 2YV pz + X{V 



v zz ) = . 



(76) 




(77) 



n=0 
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where a n and b n are real coefficients to be determined. Up to a given order M, we 
need 2(M + 1) coefficients to specify the expansion of J 7 and hence of X and Y . We 
also expand the function V about (po^o), and then substitute these expansions into 
the differential equation (176|) . For each order M in the expanded differential equation, 
there is a homogeneous algebraic equation for the coefficient of Ap q Az M ~ q , thus there 
are M + 1 separate equations. Consequently the total number of algebraic equations 
from M = up to and including order M is 

N=-{M + l)(M + 2). (78) 

However, because Eq. (175)1 is a second order differential equation, at a given order 
M, coefficients up to order M + 2 in X and Y" will appear. Therefore there will be 
P = 2(M + 3) expansion coefficients of T to be fixed at (po, z ) via this system of linear 
algebraic equations. 

For small M, N < P, and so the linear equations have solutions for any function 
V. But N grows quadratically with M, while P grows linearly, and eventually there 
will be more equations than there are coefficients of X and Y. Then, solutions for the 
coefficients of J 7 will exist only if the coefficients of the expansion of V satisfy conditions. 
The cross-over occurs when M = 4; in this case, the number of equations is 15, while 
the number of coefficients is only 14. Setting one of the equations aside for a moment, 
for the 14 x 14 system of linear equations, there will be a non-zero solution for the 
coefficients of J 77 if and only if the determinant of the matrix vanishes. The matrix 
depends only on the value of V and its derivatives up to order M + 2 = 6 at (p , z ). 
For a generic V, the determinant will not vanish, and so the only solution is for all a n 
and b n to vanish for n < M + 2. But since the expansion of the differential equation 
gives the a n and b n for n > M + 2 in terms of the values for n < M + 2, vanishing of 
the latter implies vanishing of all coefficients. Since (po, zq) was arbitrary, this implies 
that J 7 = everywhere. Even if the expansion coefficients of V manage to satisfy the 
conditions at order M = 4, they must continue to do so order by order. As a result, 
there will be a non-zero solution for J 7 if and only if V is extraordinarily fine tuned. 
The same argument can be applied to support the vanishing of the analytic functions Q 
and U in Eqs. (|52|) and (1681). 

How finely tuned must V be? In spherical symmetry, where V = V(p 2 +z 2 ) = V(££) 
and v — it is simple to show that J 7 = £ 2 is a solution to Eq. (157)1 . In this case 
the solution merely generates additional trivial Killing tensors coming from products 
involving the two additional rotational Killing vectors that exist in spherical symmetry 
(giving a total of 11 trivial Killing tensors). The Bach-Weyl solution completely lacks 
this additional symmetry. 

6. References 

[1] Will C M 2009 Phys. Rev. Lett. 102 061101 (|arXiv:0812.0110|) 

[2] Mathiina D O 2008 Integrable Systems in Celestial Mechanics (Boston: Birkhauser) 
[3] Coulson C A and Joseph A 1967 Int. J. Quantum Chem. 1 337 



Carter-like constants of motion 



12 



[4] Lynden-Bell D 2003 Mon. Not. R. Astron. Soc. 338 208 ( |arXiv:astro-ph/02104l"7j ) 

[5] Vinti J P 1959 Phys. Rev. Lett. 3 8 

[6] Weyl H 1917 Ann. Phys. (Germany) 359 117 

[7] Bach R and Wcyl H 1922 Math. Z. 13 134 

[8] Stephani H, Kramer D, MacCallum M A H, Hoenselaers C and Hcrlt E 2003 Exact Solutions of 

Einsteins Field Equations 2nd edn (Cambridge: Cambridge University Press) 
[9] Garfinkle D and Glass E N 2010 Class. Quantum Grav. 27 095004 (|arXiv:1003.0019P 
[10] Brink J 2010 Phys. Rev. D 81 022001 (|arXiv:0911.1589|l 

[11] Will C M 1993 Theory and Experiment in Gravitational Physics (Cambridge: Cambridge 
University Press) 

[12] Coclho F S and Herdeiro CAR 2009 Phys. Rev. D 80 104036 (|arXiv:0909.4413|) 



